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\ We apply Bogolubov approach to QCD to demonstrate a spontaneous genera- 

\ tion of three-gluon gauge invariant effective interaction which contributes signifi- 

cantly at the infrared region. The contribution of this interaction to a s (p 2 ) leads 
to description of its behaviour in correspondence to phenomenology and lattice 
calculations without introduction of any additional parameter. 
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In works [TJ [2] with the use of method [3] inspired by N.N. Bogolubov approach [U [5] 
a spontaneous generation of the effective interaction of the Nambu - Jona-Lasinio type 
^ ■ is demonstrated. The interaction contains no additional parameters but the QCD ones. 

Results of calculation of hadron low-energy parameters: m n , f n , m a , T a , < qq > are quite 
consistent. In the present letter we apply the approach for calculation of infrared behaviour 
of the QCD running coupling constant. 

We start with QCD Lagrangian with three light quarks (u, d and s) with number of 
colours N = 3 

3 ^ ^ 

fc=i 

C = - d u Al + g f ahc A\A c v . (1) 

where we use the standard QCD notations. In what follows we shall consider to be small 
enough and set them to zero. In accordance to the Bogolubov approach, application of which 
to such problems being described in details in work [3] , we look for a non-trivial solution of 
a compensation equation, which is formulated on the basis of the Bogolubov procedure add 
— subtract. Namely let us write down the initial expression (CQ) in the following form 



L + L 



int i 
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Lo = \ (jnM - - - A ^V^oV + ^ ■ fa bc f« v F b up f; p ; (2) 

Lint = 9s ^l,t a A^ - 1 [f; v F; v - F VFoV) - | • fabc F; v F b vp F c pfl . (3) 

Here ip is the light quarks triplet, colour summation is performed inside of each fermion 
bilinear combination, F Qpv = d p A v — d u A p , and notation ^ ■ f abc F® u F b p F pfl means corre- 
sponding non-local vertex in the momentum space 

(2tt) 4 G f abc (g pv (q p pk - p p qk) + g vp {k p pq - q p pk) + g pp {p u qk - k v pq) + 

q^Kp p - k p p v q p ) F(p, q,k)5(p + q + k) + (4) 

where F(p, q, k) is a form-factor and p, /i; q, u; k, p are respectfully incoming momenta and 
Lorentz indices of gluons and we mean also that there are present four-gluon, five-gluon and 
six-gluon vertices according to expression for F" (CQ). 

Let us consider expression (T5]) as the new free Lagrangian L , whereas expression ([3]) 
as the new interaction Lagrangian L int . Then compensation conditions (see again [3]) will 
consist in demand of full connected tree-gluon vertices, following from Lagrangian L , to be 
zero. This demand gives a non-linear equation for form-factor F. 

These equations according to terminology of works [U [5] are called compensation 
equations. In a study of these equations it is always evident the existence of a perturbative 
trivial solution (in our case G — 0), but, in general, a non-perturbative non-trivial solution 
may also exist. Just the quest of a non-trivial solution inspires the main interest in such 
problems. One can not succeed in finding an exact non-trivial solution in a realistic theory, 
therefore the goal of a study is a quest of an adequate approach, the first non-perturbative 
approximation of which describes the main features of the problem. Improvement of a 
precision of results is to be achieved by corrections to the initial first approximation. 

Thus our task is to formulate the first approximation. Here the experience acquired in 
the course of performing works [H El [3j could be helpful. Now in view of obtaining the first 
approximation we would make the following assumptions. 

1) In compensation equation we restrict ourselves by terms with loop numbers 0, 1. 

2) In expressions thus obtained we perform a procedure of linearizing, which leads to linear 
integral equations. It means that in loop terms only one vertex contains the form-factor, 
being defined above, while other vertices are considered to be point-like. In diagram form 
equation for form-factor F is presented in Fig.l. Here four- leg vertex correspond to interac- 
tion of four gluons due to our effective three-field interaction. In our approximation we take 
here point-like vertex with interaction constant proportional to g G. 

3) We integrate by angular variables of the 4-dimensional Euclidean space. The necessary 
rules are presented in papers p],[3]. 

At first let us present the expression for four-gluon vertex 

V(p, m, A; q, n, a; k, r, r, I, s, tt) = gG (f«™f*™ (u^ Z; ^ Tj 7r> A ) _ V ^ j. X , r, tt, a) - 
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—U (I, k; a, tt, t,X) + U (I, k; X, tt, r, cr) + U (p, q; tt, X, cr, t) — U (p, q; r, A, a, 7r) - 
-U(q, p; tt, a, A, r) + Z7(g, p; r, a, A, tt)) - f arn f ams (u(p, I; a, X,tt,t) - 
— Z7 (/, p; a, 7r, A, r) — U (p, Z; r, A, 7r, er) + £7 (Z, p; r, 7r, X,a) + U (k, q; tt, t, cr, A) - 
-[/ (g, k; tt, a, r, A) - U (k, q; X, r, a,n) + U (q, k; X, a, r, tt) J + 

+/ asn T mr (^(fc, P; cr, A, r, tt) - f/(p, fc; a, r, X, tt) + U(p, k; tt, r, X, a) - 
-U (k, p; tt, X, r, cr) — U (I, q; X, tt, cr, r) + U (I, q; r, tt, cr, A) — U (q, I; r, a, tt, A) + 
+U(q,l;X,a,7T,r)^y, (5) 
U (k, I; a, t, tt, A) = k a l T g wX - K l x g^ T + K h 9<rr - {kl)g UT g nX ■ 

Here a triad p, m, X etc means correspondingly momentum, colour index, Lorentz index of 
the gluon. 

Let us formulate compensation equations in this approximation. For free Lagrangian 
L full connected three-gluon vertices are to vanish. One can succeed in obtaining analytic 
solutions for the following set of momentum variables (see Fig. 1): left-hand legs have 
momenta p and — p, and a right-hand leg has zero momenta. However in our approximation 
we need form-factor F also for non-zero values of this momentum. We look for a solution 
with the following simple dependence on all three variables 

F( P i, p 2 , p 3 ) = F( Pl + P * +Pl ) ; (6) 

Really, expression ([6]) is symmetric and it return to F(x) for P3 = 0, p\ = p\ = x. We 
consider the representation ([6]) to be the first approximation and we plan to take into account 
the corresponding correction in the forthcoming study. 

Now following the rules being stated above we obtain the following equation for form- 
factor F(x) 

G 2 N ( f Y 1 f x 1 f x 

F ( x ) = / F (y)y d y - ttt- 2 / F(y)y 3 dy + — / F{ y ) y 2 d y + 

64 tt 2 \J 12 x 2 J Qx J 

x f Y „, v , x 2 f Y F(y) ,\ GgN f Y 

24 ir 2 14/16 ^(x-8y) w 9 J x/i (x-Sy) W "J 

Her x = p 2 and y = q 2 , where q is an integration momentum. We introduce here an 
effective cut-off Y, which limited an infrared region where our non-perturbative effects act 
and consider the equation at interval [0, Y] under condition 

F(Y) = 0. (8) 
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We shall solve equation ([7]) by iterations. That is we expand the last line of © in powers of 
x and take at first only constant term. Thus we have 

G 2 N ( f Y 1 f x 1 f x 

- 2 '"F (») 



y 



dy 



3G gN f Y 



(9) 



Expression ([9]) provide an equation of the type which were studied in papers [H El E] , where 
the way of obtaining solutions of equations analogous to (jUJ) are described. Following this 
way we have the unique solution of equation (Q 

i / ni at rY ^CnN C Y \ 

™ = 2 (" SWo ^ « d " + 15^" I F(V) *) (10) 

Gil ( z I? 1 /r, n _i /o _il = Co4^ z |1> 1/2) — 1) — 1/2J ; ^ = ^Q24 ^2 ■ 



'15^ 11,1/2,0,-1/2,-1 ) ~ 2z 

Constants Ci, C2 are defined by the following boundary conditions 

dF (z)- 



2 d 3 F (z) 
2z 2 + 9z 



d 2 F (z) 



2z 



dz 3 
2 d 2 F (z) 
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dF (z) 
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Here as well as in the previous papers we obtain the solution in terms of Meijer G- 
functions [6]. Conditions fl8| [TT]) defines set of parameters 

C x = 0; C 2 = 0. (12) 



OO : 



However with these parameters the first integral in diverges and we have no consistent 
solution. In view of this we consider the next approximation. We substioite solution (fIU|) 
into (jTj) and calculate the term proportional to \pz. Now we have 

F(z) = 1 + i^( lnW - 4 ''- llln(2) + ^) + 



2 r 

+ — / F(t)tdt 

3 z Jo 

4^ 



, dt 2z 



We look for solution of ({TBI in the form 

^O 2 ) — 2 ( Z 1 1,1/2, 0,-1/2,-1 



F(t) Vidt 



gVz r si( r 1 1/2 \ , 

167T 15 I 11,1/2,1/2,-1/2,-1^1 + 



(13) 



CxG^z 1 1/2, 1, -1/2, -l) +C 2 Gj°(z|l, 1/2, -1/2, -l) ; 



(14) 



Boundary conditions also change to become in explicite form 

\g%(z\1/2, 0, -1/2, 0) - (\<%(z\l/2, 0, 0, -3/2) + 

+ Gis (* I V2, o,o, 0,-3/2)) + CiGS2(*|0, 1/2, 0, -1/2,) + 
+ C 2 Gj°(z|l/2, 0, 0, -1/2,) - 

"2^ - ^(l^) + ^-4 7 -Hln(2)) =0; z = z ; (15) 
i G 30(,|l/2, 0, -1/2, l) - 9 ^ (^(,|? /2 ,o,o,i,-3/ 2 ) + 



) 2 Gil ( z 11/2,0,0,0,-3/2)) + 



+ 



15 11/2,0,0,0,-3/2 
+ C 1 Gj°(z|0, 1, 1/2, -1/2,) + C 2 Gj°(z|l/2, 1, 0, -1/2,) 

1 gVs n 

+ 47i~^r = 0; z = Zo] 

Condition F(0) — 1 leads to the following relation 

'^f**^!^ (i6 » 

which in explicite form reads 

1 + sflr 32 ^! 1 ' 1 - ^ r 32/_ ,1,1/2 \ 

\2 6 \ 2 ' 3 / 2 ' 1 ' 1 / 2,0 ' / 16 7T 26 \ '2,3/2,3/2, 1/2,0, 0^ " t " 

+ G 1 G\l(jz 1 3/2, 2, 1/2, 0, 0) + G 2 G\l(z 12,3/2,1/2,0,0)) ~~ (17) 
3 \/3^ /'l ^32/ 1 1,1/2 \ ^32/ ,1,1 



,J '9fi I " q /o 1 1 /o n n _i /o I -, ^ 1 - r 9fi * h/o 1 i n n _i /o J "I 



TT ^2 26 V ^3/2,1,1/2,0,0,-1/2; lg7r ^26^ 13/2,1, 1,0,0, -1/2, 

+ CiGi£(z 11,3/2, 0,0, -1/2) + C 2 G^(z|3/2, 3/2, 0, 0, -I/2)) =0; z = z .(18) 
We have also condition 

F(z o ) = 0; (19) 

that means smooth transition from the non-trivial solution to trivial one G — 0. The 
conditions define values of parameters 

g(zo) = 2.55779; z = 1.915838; C 1 = 0.06172743; C 2 = -0.1640803. (20) 
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We would draw attention to the fixed value of parameter z . The solution exists only for 
this value (fT2l) and it plays the role of eigenvalue. As a matter of fact from the beginning 
the existence of such eigenvalue is by no means evident. 

We consider the neglected terms of equation ([7]) as a perturbation to be taken into account 
in forthcoming studies. 

We use Schwinger-Dyson equation for gluon polarization operator to obtain a contribution 
of additional effective vertex to the running QCD coupling constant a s . The correspond- 
ing diagram is presented at Fig. 2. Due to this vertex being gauge invariant, there is no 
contribution of ghost fields. So the contribution under discussion reads 

An ( , 9 gn f v%M -g - f , g - 1) rjttth i±izi±f F ^ q2 ± £ dq ^ 



2(2 7r) 4 y (g 2 + p 2 /4) 2 - (pq) 

Using expression (j3J) after angular integrations we have 

AU^(x) = (g^ v p 2 - VyPv) n(ar) ; x = p 2 ; y' = q 2 + — ; 

n( ,) = .!™f'f ( 16 ^ _ 48 ^ + «» - | + (22) 

v; 32 vr 2 Vx 2 y 3x/ 4 y'-x/2 V x 2 x y 2 1 



Y 



Fjy'W 
L y'-x/2 



Here coupling (7 corresponds to Integrals of Meijer functions depending on y' multiplied 

by any power of y' can be evaluated analytically. In view of this we expand the denominator 
in (|22p in series in powers of x/2. Substituting solution ffTU]) into integrals we obtain explicit 
expressions for each term of expansion of expression f[2"2"j) . Calculations show that results 
with different numbers of expansion terms begin practically coincide starting from three-term 
expansion. In what follows we present results with four terms. 
So we have modified one-loop expression for a s (p 2 ) 

a,(x) = : , 9N : ; ,/ , z— - ; X = p 2 J 60 = 11 L . 23 

v; 47r + 6oa 8 (pg) ln(x/A 2 ) + 4vrn(x) ' y% 3 v 7 

It is remarkable that function Il(x) at x = Y (z — z ) turns to be very small, almost 
zero. We just expect this quantity to be zero exactly. So this property of the approximated 
polarization operator indicates the consistency of the procedure being used. Namely we 
normalize a s {p 2 ) at point po, which correspond to our cut-off Y. Coupling constant g entering 
in expressions f)16p and (I2"2"j) is just corresponding to this normalization point. Performing 
the well-known transformations in expression (I23p we have for u < Uq 



4tt / , s 8 VNn \ x Y 



a s (u) = —(\n(u) + — 1 J ; u = — 2 ; u Q = — 2 = 14.61328; 

bo \ b g J A 2 QCD A 2 QCD 
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I = - 



F (t)dt 



16- 



£ 3/2 



t 
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F (t)dt 



48^= + 45Vt 

>z 2 

G 2 Nx 2 



Vi(2Vi - yfz) _ 
For u > Uq we use the perturbative one-loop expression 

47T 

Oo ln(it) 



1024 7r 2 



+ 



&o = 9. 



(24) 



(25) 



Two expressions (1251) . ( 1241 give equal values ao for a s at point zq, while value a s (uo) (1201 is 
defined by equation ( |T6l) . The self-consistent result for expressions ( 1241 [251) with account of 
previous relations is unique and reads as follows 



2 



1.915838 



"o 



d = 0.06172743; C 2 



14.61328; a = 0.5206188 
= -0.1640803. 



(26) 



Behavior of a s ( 1241) with u = Q/Aqcd is presented at Fig. 3 for Aqcd = 0.2GeV and 
0.05 GeV < Q < lGeV. The behaviour with maximum at Q ~ 0.6 GeV and maximal 
value a^ a:r ~ 0.55 agrees to calculations in work [7]. Qualitately the result also corresponds 
to lattice calculations in work [8] (see also discussion in paper [9]). Note, that we begin plot 
at Fig. 3 starting from p = 0.05 GeV, because a s (u) has a pole at very small u, which is 
analogous to the well-known perturbative pole at u = 1. Now this pole is shifted to the far 
infrared region. One may deal with it using the method proposed in work [TU] and subtract 
from (1241) the following term 



47T 



b D (u - mqo) 



«00 



0.005769; D = 170.1594. 



(27) 



This procedure practically does not change the result presented at Fig. 3 in the denoted 
interval of Q. For comparison we present the modified a s {Q) for interval 0.01 GeV < Q < 
1 GeV at Fig 4. Value of a s (Q) at zero reads a a (0) = 1.4205. 

All values are now expressed in terms of Aqcd- Emphasize, that there is no additional 
parameters to describe the non-perturbative infrared region. 

From definition of variables z, u and x we obtain value of new coupling constant G 



G 



5.497571 

1V QCD 



(28) 



We calculate non-perturbative vacuum average of the third power in gluon field, which 
is immediately connected with our results. We have 



< 9' 



f jp a TP" ip c 
J abc r r vp r pfi 



> 



g 3 G96n A 
(2tt) 8 



hh 



110.77 A' 



QCD i 



(29) 
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h 



I 




! (i _ y )2(i _ 32/ /4) 



0.278756 ; I 2 



I 



Y 



x 3 dxF(x) = 0.18062 



y(i-y/2) 

(32 tt) 4 
2G 4 iV 2 ' 



Value ( 1291) seems reasonable. However the accuracy of phenomenological and lattice defini- 
tions of this value is not yet satisfactory. 

To conclude we would state, that method [3] being applied to non-perturbative a s proves 
its efficiency even in the first approximation, which is considered here. Bearing in mind also 
results of works [TJ [2] on application of the approach to low-energy hadron physics we would 
express a hope, that in this way we could obtain the adequate tool to deal non-perturbative 
effects in QCD and, maybe, in other problems. 

The mean non-perturbative value for a s turns to be here around 0.5. For this value 
results of works [H[2] seem to be quite consistent. 

The author express his gratitude to D.V. Shirkov for valuable discussions. 
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Figure captions 



Fig.l. Diagram representation of the compensation equation. Black spot corresponds to 
tree-gluon vertex with a form-factor. Simple point corresponds to a point-like vertex. In- 
coming momenta are denoted by the corresponding external lines. 

Fig. 2. Loop contribution to polarization operator. 

Fig.3. Behaviour of a s (Q), for 0.05 GeV < Q < 1 GeV, A QCD = 0.2GeV. 

Fig.4. Behaviour of modified a s (Q), for 0.01 GeV < Q < 1 GeV, A QCD = 0.2GeV. 
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Fig. 2. 
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